Memory function approach to the Hall constant in strongly correlated electron 

systems. Part II 
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The anomalous frequency and doping dependence of the Hall constant in the normal state of 
high-Tc superconductors are investigated within models of strongly correlated electron systems. In 
Mori theory, the transition of the Hall constant from infinite to zero frequency is described by a 
memory function. It naturally introduces a second time scale, that, within the t — J model, is 
identified with the spinon relaxation time of Anderson. This provides us with a phenomenological 
understanding of the interplay between the frequency and temperature dependence of the Hall 
constant for frequencies below the Mott-Hubbard gap. Within the single-band Hubbard model in 
the limit U 2> t (U: Coulomb repulsion, t: hopping amplitude), the memory function is calculated 
via its moments and shown to project out the high-energy scale U, introduced by doubly occupied 
sites. This causes the Hall constant to decrease by a factor (l + <$)/2 (S: doping), when the frequency 
is lowered from infinity to values within the Mott-Hubbard gap. Finally, it is outlined, how the 
Hall constant may be calculated in the low frequency regime. 



I. INTRODUCTION 



Recently, we have proposed a memory function treat- 
ment of the Hall effect in strongly correlated electron 
systems [Q. Its main advantage consists in the fact that 
we are directly provided with a representation of the Hall 
resistivity. Therefore, in contrast to ordinary approaches, 
we may dodge the issue of coping with a quotient of con- 
ductivities, which may be precarious due to resonances 
like the Drude peak. 

In our approach, the frequency dependent Hall con- 
stant is given as the sum of its infinite frequency limit 
and a memory function contribution. A perturbative cal- 
culation of both terms within the single-band Hubbard 
model has demonstrated the usefulness of this approach 
[jlj . Our results turned out to be in qualitative agreement 
with the unusual experimental findings for high-T c super- 
conductors as, e.g., La2-iSr,5Cu04, except for the doping 
dependence of the Hall constant in the vicinity of half fill- 
ing, i.e. below S ~ 0.1. There, a 1/5 law is observed which 
requires a description in the strong-correlation regime Q . 
This parameter regime will be the subject of the present 
paper, which is organized as follows. 

In Sec. II, the main results concerning the representa- 
tion of the Hall constant within Mori theory are compiled 
from Ref. ffl. Then, in Sec. Ill, we study the moments of 
the memory function in the limit U — > oo. Thereby, the 
memory function is found to eliminate the high-energy 
scale set by the Hubbard repulsion U. This provides us 
with an explanation of the frequency dependence of the 
Hall constant in the cross-over regime from w > [/ to 
W <C u) <C U, where W is the bare band width. Fur- 
thermore, a simple analytical treatment of the Hall ef- 
fect within the Hubbard I approximation is presented, 
the results of which turn out to be in full accord with 
that obtained by high-temperature expansions. In Sec. 
IV, we reformulate the frequency dependent Hall con- 



stant within the t — J model in order to address the low 
frequency regime as well. The comparison to a repre- 
sentation of Rh(lo) in terms of two relaxation rates and 
effective masses introduced by Anderson || provides us 
with an interpretation of the additional relaxation rate 
th in the language of the Mori theory and a phenomeno- 
logical understanding of the anomalous frequency and 
temperature dependence of the Hall constant in high-T c 
superconductors. Finally, we explain how the emerging 
picture may be put onto the basis of a microscopic cal- 
culation. In this context, we reduce the problem of cal- 
culating the memory function in the t — J model to the 
much easier one of finding the first few moments of the 
ordinary current-current correlation functions. We con- 
clude this paper in Sec. V with a discussion of our main 
results. 



II. MEMORY FUNCTION APPROACH TO THE 
HALL CONSTANT 

For a more detailed introduction into the following 
memory function treatment of the Hall effect in strongly 
correlated electron systems, the reader is referred to Ref. 

0- 

We choose the simplest model of strongly correlated 
electrons, namely the single-band Hubbard model on a 
two-dimensional square lattice: 



H 



E 

<ij>cr 



Pij C i<T C 3 



(1) 



The second term is a local repulsive interaction. The sum 
in the hopping term is restricted to nearest neighbors and 
the Peierls phase factor ~ exp{ieA(Rj)(Ri — Rj)} 
guarantees the gauge invariance Q (the sign of the charge 
e is chosen to be negative). The vector potential decom- 
poses into two terms describing the electric and mag- 
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netic field, respectively: l(t,r) = A el {t) + A mag (f), 
E(t) = -jj-A el (t), and H = rotA ma s(r). In linear re- 
sponse theory, the current operators are defined as: 



Ju ■ = 



1 5H(t) 



(2) 



A eJ =0. 



A suitable representation of the Hall constant is obtained 
within Mori theory. This formalism enables us to sepa- 
rate the dynamics of the current operators J x and J y from 
that of all the other degrees of freedom, the influence of 
which is then accumulated in memory functions 0. In 
the center of this separation stand the superprojectors 



= i-p. 



(3) 
(4) 



acting in operator space C and being defined with respect 
to Mori's scalar product 



two time scales with the relaxation rates of the holon 
and spinon degrees of freedom in Anderson's tomographic 
Luttinger liquid theory ||] . The memory function (||) can 
be represented as a spectral integral 



M(z) = J *± 



dto M"(lo) 



(10) 



where the spectral function M"(oj) is given by the dis- 
continuity across the real axis: 



M(lo ± i0+) = M'{oj) ± iM"(uj) , 



(11) 



and can be shown to be real and even. Since the mem- 
ory function vanishes as 1/ z 2 at high frequencies, the first 
term on the Rhs. of Eq. ([?]) represents the Hall constant 
in the limit uj — > oo, that was considered by Shastry et 
al. [0. It will be denoted as from now on. 



III. HALL CONSTANT IN THE CROSS-OVER 
REGIME FROM tj > £/ TO to < U 



(A\B) := (1/(3) / dr < exp{rL}A + B > . (5) 
Jo 

Here, (3 is the inverse temperature, the Liouville opera- 
tor L maps a given operator onto its commutator with 
the Hamiltonian: LA = [H,A], and x° — P(Jx\Jx) is the 
static current susceptibility. This scalar product implies 
the so-called Kubo identity 



P(A\L\B) =< [A+,B] > 



(6) 



which will play an important role in the following sec- 
tions. By representing the current-current correlation 
function within the Kubo formula for the conductivity 
tensor in terms of a memory matrix, we obtain the fol- 
lowing representation of the Hall constant: 



Rh(z) 



N CI 
hm — 



iM(z) 



ie 2 x° 



H 



(7) 



The first term is given by the cyclotron frequency CI — 
(l/x°) < [Jx j Jy] > while the second is a memory func- 
tion contribution: 

M(z) = £(QLJ x \—^—\QLJ y ) . (8) 

Here, z is a complex frequency that ultimately has to be 
continued analytically to z = u> + i0 + . This function has 
the structure of a relaxation function for the so-called 
residual forces 



fv — iQLJu 



(9) 



whose dynamics is governed by the projected Liouville 
operator QLQ rather than L Thus, these forces may 
vary on a time scale that is differ ent from that of the cur- 
rent operators J v . In Sec. IV A, we shall identify these 



Numerical studies of multi-band Hubbard models for 
the Cu-0 planes of high-T c materials indicate that all 
low-energy excitations are reproducible within a single- 
band Hubbard model with U ~ W (]]]. Unfortunately, 
this parameter regime is not accessible to reliable ana- 
lytical calculations. Therefore, we shall exaggerate the 
impact of U by considering the range U 3> W instead. 



A. The moments of the memory function in the 
limit U — > <x 

The main problem in dealing with the memory func- 
tion (||) is related to the fact that its dynamics is governed 
by the projected Liouville operator QLQ rather than L. 
We resolve this difficulty by inquiring into the properties 
of the memory function via its moments. Then, we may 
get rid of the superprojectors Q simply by resorting to 
their definition through Eqs. (||) and (^). 

We start by writing the Hall constant (|t|) in terms of 
spectral functions: 



Rh{z) 



N 



e 2 ( X °) 2 H y_ 



duj \ S h (lu)+K{lu) ■ 



(12) 



The first spectral function in the integrand is defined as 

S B (w) = -i<[js,8(u + L)J v ] >, (13) 

and corresponds to the Hall matrix element of the 
current-current correlation function. The second func- 
tion is related to the spectral function M"(ui), introduced 
in Eq. (|l]), by K(uj) = ^M"(u})/(wu), which implies: 

K{u) = —(QLJ x \5(uj + QLQ)\QLJ y ) . (14) 
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This function describes the correlation between the resid- 
ual forces f x and /„. In the following, we shall unravel a 
simple connection between the moments of the two func- 
tions Sh{u) and K(u>) in the limit U — > oo. 

First, we note that both functions may be shown to be 
even and real fl]J|. Furthermore, we expect both func- 
tions to vanish beyond a certain frequency above U, since 
U is the highest energy scale in the problem. This as- 
sumption will be corroborated below up to mistakes of 
the order t/U. Hence, for finite but large U, all moments 
exist and it is sufficient, to consider only the even ones: 



/oo 
dioS H {uj)oj 2 \ 
-OO 
/"CO 

K, = I <LjK{oj)uj 2 \ 



(15) 
(16) 



where I > 0. Since we are mainly interested in the dc- 
Hall constant, we would like to calculate K$. Unfortu- 
nately, this is not feasible on the basis of Eq. ( |l4| ) , since 
the inverse of the projected Liouville operator QLQ does 
not exist. However, all the other moments ( |l6| ) can be 
calculated: Due to Eq. (|l4|), they are given as: 



Ki = i(3(J x \ LQL . . . LQL \J y ) (Z > 1) 
21 projectors Q 



(17) 



Now, if we insert Q = 1 — P, this expression decomposes 
into 4 Z terms. Consider a special one consisting of p su- 
perprojectors P. By using the definition of P, this term 
is seen to decompose further into / = p + 1 factors of the 
form P(J\L mj \J), the orders of which are U mi ~ Y due to 

the Kubo identity (g). Here, J2j=i m j = 2Z + 1. Thus, 
each factor lowers the relative order in U by one. Hence, 
the more superprojectors P a given term is composed of, 
the lower its relative order in U is. Therefore, all super- 
projectors Q may be removed from the Rhs. of Eq. ( |l7| ) 
to leading order in t/U. And this, in turn, establishes 
the following relation between the I > 1 moments of Eqs. 
tfTJ) and ©: 



jj2i u 21 [ U h 



(18) 



This equation does not imply that the functions (|T^) and 
differ only by a sign in the limit U — > oo. This con- 
clusion would require positive or negative definite func- 
tions and finite moments even in the limit U — ► oo. 
None of both conditions is satisfied. To proceed any- 
way, we remind ourselves that in the context of the Hub- 
bard model in the strong correlation limit, any spectral 
function is believed to separate into individual "peaks" 
centered around integer multiples of U ||-[l(| (in this 
context, any connected structure of a given spectral func- 
tion, irrespective of its detailed shape, is referred to as 
a "peak"; for instance, it may vanish at discrete points). 
This reflects the fact that one-particle excitations may 
be grouped into two Hubbard bands separated by the 



so-called charge transfer gap, which is of the order U. 
Since the current operators produce particle-hole excita- 
tions, we expect the functions (|l3|) and ( |l4| ) to have peak 
structures centered around uj — and ±U related to ex- 
citations within the two Hubbard bands and across the 
charge transfer gap, respectively. Therefore, these peaks 
are expected to have widths of the order of those of the 
Hubbard bands. In the following, we shall prove this pic- 
ture at least for the function ([l3|) and derive formulas 
required to extract more information from the relations 



B. Structure of the functions Sh(w) and K(ui) 

The appropriate technique to investigate spectral prop- 
erties of the Hubbard model in the strong correlation 
limit was pioneered by Harris and Lange in the special 
case of single-particle excitations || and generalized by 
several other authors, see for instance Refs. Jj|[l0|]. At 
the heart of this procedure stands the decomposition of 
a given operator into terms, which increase the number 
of doubly occupied sites by integer values p: 



O 



(19) 



Together with the Lehmann representation of a given 
spectral function, one may then address the proper- 
ties of its individual peaks. The decomposition jl9| ) 
is accomplished by an iterative procedure based on a 
canonical transformation of the Hubbard Hamiltonian: 
H — > exp{iS}H exp{— iS}. The expansion of the opera- 
tor S up to the order I in t/U eliminates those processes 
from H which change the total number of doubly occu- 
pied sites up to the order t l /U 1-1 . The corresponding 
transformed Hamiltonian H^ l+1 "> in turn helps to fix the 
next order of S in t/U and so on. Thus, subsequent 
iterations generate increasing orders in t/U. Once the 
generator S has been found to a given order, one may 
decompose any operator to the same order by first de- 
composing its rotated counterpart [ ^0[ . 

In our case, t/Z7-expansions may be terminated after 
the zeroth order term since the relation (18) indicates 
that it is not sensible to go beyond. Then, we do not have 
to distinguish between original and transformed Fermi 
operators and the decomposition (|l^) specialized to the 
case of the operators = cf a Cj a , making up the com- 
ponents of the current operator, becomes: 



D" 



D 



ir-U 



D 



ij-0 



D 



In terms of Hubbard operators X, 



ij;U 



Oct 



(20) 



(1 



(1 - h is )cj a , X 



n-2 



' and X? s 



' H(T } 1 



X, 

where a = —a, the terms of Eq. (|20|) may be written 
conveniently as: 



3 



U ij-U 



X a ® X^ a -4- x 2a x a2 



2a v 0a 



xrx 



(21) 
(22) 
(23) 



The Lehmann representation of the function ([L3|) is de- 
rived straightforwardly: 



S ff (w) = !^M nm (e-^" 



) 



x5(w - [e n - e m ]) , 

{< n\j x \m >< m\j y \n > 



2i 



< n\J y \m >< m\J x \n >} . 



(24) 



(25) 



where states and energies are defined through the eigen- 
value equation (H — /iN)\n >= e n \n >. Inserting the 
decomposition of the current operators corresponding to 
( pp| ) , we find that the peak centered around u = pU has 
a weight given by 



P u 



2i 



{< [Jx-pUi Jy-plj] > + < [Jx;- 



■pU i Jy;pU . 



>) 

(26) 



and that only the peaks p — and p = ±1 survive in 
leading order in t/U . 

We assume that the function has qualitatively the 
same triple-peak structure. Although not proven, this 
assumption is shown to lead to reasonable conclusions. 



C. Frequency dependence of the Hall constant in the 
range uj 3> W 

Given the peak structure of the functions (|l3| ) and 
( |l4| ) , only the contributions of the satellite peaks around 
uj — ±U can be resolved in the I > 1 moments (|l5|) and 
(|l6|). Thus, the relation ( |l8| ) implies that the "spectral 
weights" of the peaks of K {uj) and Sh(u>) around uj = U 
differ only by a sign. Together with Eq. (fL2|) , we may then 
draw the following conclusions: For u ^> U, all peaks of 
Sh(uj) contribute to the Hall constant and none of K(uj). 
In the frequency range W <C uj <^ U, the high-frequency 
peaks cancel each other out while the contribution of the 
zero-frequency peak of K(uj) is negligible. Within this 
charge transfer gap region, the frequency dependent Hall 
constant is then lowered by a factor p in comparison to 
its infinite frequency limit, 



R* H = R H (W < u) < U) = pRS, 



(27) 



if we define p to be the relative spectral weight of the low 
energy structure of Sh(uj): 



Wn 



_ - [^r;0, Jy,o] 

W-u + W + W V ~ < [J x J y ] > 



(28) 



Eq. ( p7| ) is valid except for mistakes of the order t/U. 
Therefore, it is sufficient to evaluate it in the limit 
U -> oo. 

First, we seek for an interpretation of R* H - We begin 
by noting that the infinite frequency Hall constant, i.e. 
the first term of Eq. (0), may be rewritten as follows: 



1 2a d 



e (2a„) 



(29) 



where a n and is the amplitude of a nearest neighbor 
hop and a hop diagonally across the unit cell, respec- 
tively. This holds on account of the following relations, 
which may be proven by straightforward analysis: 



< [Jx,Jy] > = SiNet 2 Ha d , 
X° = AtNa n . 

In the limit U — > oo, we have 



lim a d = < XfX"^^ >. 

U^oo R R+x+y 



lim a n = <XfX% +& > 



(30) 
(31) 

(32) 
(33) 



where x and y is a primitive lattice vector in the x- and 
y-direction, respectively. Here and in the following, ex- 
pectation values are taken with respect to states without 
double occupancies. Then, the projected current opera- 
tors appearing in the nominator of the Rhs. of Eq. fl2S| ) 
take on the following form: 



J v . = it ^2 



Av r> vcrQ vOcr 



(34) 



<ij>a 



since the term Xf^XJ 2 of Eq. ( |2l| ) may be omitted. Here, 

Aij = Ri — Rj and Py is the phase factor defined in the 
text following Eq. (|l|). From Eqs. (pq)-(|3^), we conclude 
that R* H , defined in Eq. (CT), represents the infinite fre- 
quency Hall constant of the U = oo— Hubbard model, 
which is defined as follows: 



H 



-t 



E 

<ij>a 



p. l3 xfxf . 



(35) 



In fact, the analysis of Sec. II is straightforwardly car- 
ried over to this model, the current operator of which is 
then found to be given by Eq. Note, that p ^ 1 

expresses the fact that the limits u> — > oo and U — > oo 
do not commute: If we start with the limit uj — > oo, the 
integral over Sh(lo) in Eq. ( |T2| ) extends over all three 
peaks while when taking the limits in reversed sequence, 
the high-energy peaks are unattainable from the outset. 

Next, we derive an exact analytical expression for p 
by taking the additional limit T — > oo. We may won- 
der whether this is reasonable. However, since the limit 
U — > oo was already carried out, at least, the condition 
T <C U is satisfied, i.e., the thermal energy cannot excite 
an electron across the charge transfer gap. Furthermore, 
we expect neither of the high-frequency objects R^ and 
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R* H to depend appreciably on temperature, since they 
correspond to and generalize the semi-classical expres- 
sion for the Hall constant In the context of high- 
temperature expansions, one has to cope with electrons 
or holes hopping around closed loops, which are defined 
by a sequence of adjacent lattice sites. Therefore, it is 
convenient not to expand the phase factors. Then, an 
electron hopping along a polygon ijk . . .li accumulates 
a phase proportional to the flux 4>ijk...u enclosed. The 
procedure to expand expectation values of Hubbard op- 
erators like that of Eqs. ( |32| ) and ( p3| ) or the nominator 
of Eq. ( |28| ) in powers of 1/T is explained, e.g., in Ref. 
prf. To leading order, we obtain: 



P 



1 + 5 



(36) 



where here and in the following, the electron density is 
measured via the average number 5 of holes per lattice 
site introduced into the half filled system, i.e. 6 = cor- 
responds to half filling while 6=1 corresponds to the 
empty band. Although the high-temperature calculation 
has introduced another high-energy scale into our sys- 
tem, we expect this result to hold qualitatively for low 
temperatures as well. For instance, since p is a measure 
for the difference between the plateau values of Rh{oj) 
on both sides of U , it is expected to increase monotoni- 
cally as half filling is approached. Further down, we shall 
derive the same expression for p within a simple approx- 
imation valid at T = 0. 

Finally, we calculate to leading order in 1/T. The 
leading order of the amplitudes ( |32| ) and (|3^) are found to 
be a n = (pt/2)6(l-6) and a d = -{(3 2 t 2 /4)5{1-5){1-Z5), 
which results in 



Id 



1 



1 



(37) 



Together with Eqs. ( f27| ) and (pq), we recover the result 
for R* H of Ref. |], that was derived within the t — J 
model in leading order in 1/T. In this work, it was fur- 
ther shown that, although R* H is renormalizcd as a func- 
tion of T and J when including higher orders in 1/T, 
the doping dependence of the Hall constant retains its 
most important features: its sign change at 5 ~ 1/3 and 
its singular behavior in the vicinity of half filling. What 
seems to be striking at first sight is the fact that the same 
properties are encountered for Rff, i.e. the Hall constant 
at frequencies well beyond U. At such high frequencies, 
the dynamics of the electrons is insensitive to the interac- 
tion U. However, for nondynamical quantities as matrix 
elements, the correlations introduced by the Hubbard in- 
teraction remain important. 

Before we discuss possibilities to extend the moments 
technique to lower frequencies, we shall rederive the re- 
sults of this subsection within a simple approximation, 
valid at T = 0. 



D. Hubbard I approximation 

An expression for the frequency dependent Hall con- 
ductivity with vertex corrections having been neglected 
was derived in Ref. @: 



e 3 H 



E 



ka 



dk r 



2 d\ U H (z,k) 
~dkj ' (38) 



<?xy(z) 

n ff (iw m , £) = 4 E G k( iuj n) 2 

n 

x [G^{iuj n + iu> m ) - Gg(iu n - iu> m )] . (39) 

Here, the Green's function is given in terms of its spectral 
function as 



duj AAw) 



2tt iui n 



(40) 



and e, in our notation, is negative. Furthermore, it is as- 
sumed that the momentum dependence arises solely from 
the dispersion eg of the bare band. In Mori theory, the 
Hall conductivity is represented as M: 



<J X y(z) 



w- 



L 



Jy). 



(41) 



Therefore, the function ( |l3| ) may be connected to the 
Hall conductivity <t xv (oj ± i0 + ) = o' h {oj) ± ia'^(uj) via 



S h (oj) = 



(42) 



Since cr'jj(u)) arising from Eqs. ( |38| ) an d (|39| ) may be 
shown to be real and odd, the function (|13|) has indeed 
the correct analytic properties. Inserting Eqs. ( |39| ) and 
@ into Eq. @ and using Eq. (f42|), we obtain, after a 
standard calculation [|13l : 



dk x 



2 «2 



(43) 



Ag(e - w) 
2tt 



UJl — U! 2 

[/(e) -/(e-u,)] 



+ (uj — > —w) 



(44) 
(45) 



By means of a partial integration, the corresponding sum 
rule is straightforwardly shown to be satisfied (cf. Eqs. 
©, ©, and ©). 

Next, we evaluate Eqs. (^)-(^) in the so-called Hub- 
bard I approximation 111]: 



Ar(u) 1 + S 



8{uj- 



1 



2tt 2 

l-i!., TT 1-5 . 
+ 6(w-U- -^-cg). 



(46) 
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At U — oo, only the first term contributes. Then, the 
quantity ( ]44|) becomes 



XAu) = 2 



(47) 



and by means of an integration by parts, we obtain 



Sh(oj; U = oo) 



l + S 



-i<[J x ,J y ]>)6(w). (48) 



Here, the expectation value is given by Eqs. (|30|) and 
(pl|). The sum rule obeyed by Sh{^) for U ^ oo im- 
plies p = (1 + S)/2, as in Eq.[36]). This provides fur- 
ther evidence that the result (|36) may be trusted for all 
temperatures. Note, that in a real system, the function 
Sh(<-o) must vanish at lo = 0. Only then do we obtain 
a finite Hall conductivity at u> = 0. The result (^|) is 
an artifact of the vanishing width of the lower Hubbard 
band in (46). If we consider finite but large values of U, 
we may simulate the high-energy peaks of the function 
(O) by delta functions as well: 



Sh(w, U — > oo 
'1 + 8 



-l < [J X , Jy] > 

1-8 



[6(w -U) + 5{u> + U)]\ . (49) 



This may be proven by calculating the moments 
the limit U — > oo. 

Finally, we calculate R^, given by Eq. (^), analyti- 
cally within the approximation ( |46| ) and at U = oo. Al- 
though this has been done numerically some time ago 
p5[ , our simple analytical treatment allows for a direct 
comparison with the exact high-temperature result (p^) 
and demonstrates that the resulting doping dependence 
does not rely on the location of the Fermi surface. This 
last point is blurred in the Boltzmann equation based 
approach of Ref. Jl6| . 



The amplitudes on the Rhs. of Eq. (29) may be written 
as follows: 



a,, 



(50) 
(51) 



Here, the density is given as 

1 + 5 „,l + 6 



nr = 

2 



(52) 



since in Hubbard I approximation at U — oo, only the 
first term on the Rhs. of Eq. ( f46| ) is to be kept. Momen- 
tum dependences arise solely from the dispersion of 
the bare band, why next, we are looking for expressions 
for the following functions: 



fc 

B(e) = — ^2 cos k x cos k y S(e- eg) 



(53) 
(54) 



A convenient way to find smooth approximations for 
these functions is to calculate them in the limit of in- 
finite dimensions d 101: 



A(e) = -^=D°(e), 
Bte) = \^-\)D\e), 



(55) 
(56) 



where D°(e) — (1/y/ir) exp(— e 2 ) is the density of states 
of the bare system. Except for the prefactors 1/y/d and 
1/d, which ultimately cancel each other out, the func- 
tions (p5|) and (|6|) may be compared to those calculated 
numerically on a two dimensional lattice. This reveals 
that the main effect of the limit d — > oo is to smooth out 
the logarithmic singularities at zero ener gy e ncountered 
in the case of the functions D°(e) and (|54| ) in d = 2. 
Therefore, this limiting procedure does certainly not af- 
fect the validity of our present analysis in any serious 
manner. At T = 0, Eqs. ©-© imply: 



2a T 



2a d 



1 



2Vd 
-e F 



y i + s' 



(57) 
(58) 



A relation between the Fermi energy and the doping pa- 
rameter is established straightforwardly, which, in terms 
of the function 



H(S) 



2e F (6) 
l + S 



may be written as: 
1 - 



l 



erf (7J(c$)) 



(59) 



(60) 



In summary, the function 
be written as follows: 



9h at U = oo and T = may 



H(S) 



l+S D°(H(5)) 



(61) 



Firstly, at S = 1/3, 
we recover the exact 



The most important features are: 
Rfl vanishes. Secondly, at 5 — > 1 
result i?f ~ — l/(je| (1 - S)). And thirdly, at S -► 0, we 
find ~ 1/(| e| 25). The last two statements are proven 
with the asymptotic relation D°(e) /e ~ ±1 — erf(e), valid 
in the limit e — > ±oo. All these points are in exact agree- 
ment with (^7|). We take this as an indication that, on 
the one hand, the high-temperature result ( |37| ) remains 
qualitatively valid even at low temperatures, and, on the 
other, that the Hubbard I approximation is remarkably 
good in the case of the quantity (p9|). In Ref. p6l, the 
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doping dependence of the Hall constant in Hubbard I ap- 
proximation was discussed in terms of the Fermi surface. 
This is misleading for two reasons. 

For one thing, the Hubbard I approximation misplaces 
the Fermi surface: The Luttinger theorem, which relates 
the volume enclosed by the Fermi surface to the electron 
density fl7j| , is violated in this approximation. In con- 
trast to this, angle resolved photoemission experiments 
(ARPES) on cuprates like Ndz-zCe^CuC^ gl,@ ap- 
pear to be consistent with LDA bandstructure calcula- 
tions which, in turn, imply the validity of this theorem. 
Despite this flaw, the approximation (|46| ) yields a dop- 
ing dependence for the Hall constant which is in good 
agreement with experiments on I^-zSr^CuCj. 

For another, it was pointed out in Refs. || and [po[ , 
that the high-frequency object is not directly related 
to the location and topology of the Fermi surface. In- 
stead, in a strongly correlated system, the entire Bril- 
louin zone tends to get populated. In consequence, the 
weighted density average ( pl|) receives contributions from 
the entire Brillouin zone rather than just from the vicin- 
ity of the Fermi surface. 

We can demonstrate this more explicitly by slightly 
changing the form of the lower Hubbard band in Eq. 
(fl6|) : We broaden the delta function a little bit and shift 
some relatively small amount Z of spectral weight to a 
new delta function contribution ZS(lu — L(eg)), with the 
function £(eg) being undetermined yet. This amounts to 
replacing Eq. (|52| ) by = h{e^) with the function h(ej;) 
differing from ^rf{^^- e £) only in the following respects: 
The step at eg = |S an d of height -Mp is smoothed out 
while a new, much smaller one of height Z occurs at 
eg = L _1 (ei?)- This last condition fixes the new Fermi 
surface. By choosing the function L(e) appropriately, we 
may place the Fermi surface wherever we want. As long 
as Z < irp, the crucial average ad — dcB{e)h{e) 



does not differ very much from the result in Eq. (58). 
This reasoning illustrates that, in the presence of strong 
correlations, the doping dependence of is in fact un- 
corrected to the Fermi surface location. Also note, that 
it does not matter whether the function h(e) arises from 
coherent or incoherent excitations. 



IV. HALL CONSTANT IN THE LOW 
FREQUENCY REGIME 

In this section, we discuss the frequency dependence of 
the Hall constant for frequencies below the Mott-Hubbard 
gap. Therefore, the appropriate model to start with is 
the t — J model. It is straightforward to show that Eqs. 
(p^|)-(p^) are still valid, however, with all quantities be- 
ing redefined within the t — J model . Apart from the 
redefinition of the Liouville operator, this amounts to re- 
placing the canonical Fermi operators through projected 
ones in all quantities that appear, i.e. ci a — ► Xf a , and 
c+ — > X?°. In particular, the current operator is then 



given by Eq. (|34|). If we renormalize the functions ( |13| ) 
and (jlj) according to 

S H (u; t — J) = —i < [J x;0 , Jy. ] > s(oj) , (62) 
K(ui; t- J) = -i< [J x . , Jy-o] > k(u>) , (63) 
the analog of Eq. (113) reads: 



R H (z)=R* H 1 



dtv k{u>) ■ 



(64) 



Here, R* H is the infinite frequency Hall constant of the 
t — J model that was already investigated in Ref. ||, 
and which has been introduced in Eq. ( p7j ) in the special 
case J — 0. Furthermore, we have taken into account 
that the function s(u>) is normalized to unity, while k(u>) 
represents the unknown me mory function contribution. 



HI Bj , we expect the func- 



From the discussion in Sec. 

tion k(u>) to have only one peak centered around zero 
frequency, because, in the t — J model, doubly occupied 
sites can occur only virtually. 

Before we set about discussing possibilities to calcu- 
late this function via its moments, we try to gain some 
phenomenological insight. 



A. Phenomenological discussion 

Very recently, the normal state ac-Hall constant was 
measured in YBa2Cu3C>7 thin films for frequencies up 
to 200 cm" 1 p2[ . In this work, the experimental data 
have been fitted successfully in terms of parameters in- 
troduced by Anderson || to account for the observed T 2 
dependence of the inverse Hall angle in high-T c mate- 
rials. Anderson's theory is based on spin charge sepa- 
ration with two different relaxation times and effective 
masses associated with the spinon and holon degrees of 
freedom: r tr is the decay time of the holons with effective 
mass mtr, scattering off thermally excited spinons. On 
the other hand, a transverse relaxation rate 1/th is de- 
termined by the scattering between the spinons. Apart 
from this, o~ xx and o~ xy have the ordinary Drude form, i.e. 
o~ X x oc Ttr/mtr and a xy = g xx uj c th- Here, the cyclotron 
motion is characterized by a mass i.e. oj c cx 1/m^f. 
In Ref. [p2| , Anderson's theory was extended to finite fre- 
quencies via the replacements r tr — > Tt r /(1 — iu>T tr ) and 
tji — > tjj I (1 — iojTii)- This led to the following represen- 
tation of the frequency dependent Hall constant: 



R H {(J) = 1 + : 

inn ne \ T tr 1 — iujth 



(65) 



This result is equivalent to the exact expression (|64|), 
provided the following identifications are made: 



k(uj) = — — — L 1/th {lj) , 

~tr 

D* m tr 1 
% = ■ 

mn ne 



(66) 
(67) 
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Here, Lr(u>) denotes the Lorentzian of width T normal- 
ized to unity. Therefore, the unusual relaxation time th 
is a measure for the width of the function (^33|) and thus 
determines the decay rate of the correlation between the 
residual forces f x; o and f y -o of the t — J model (cf. Eqs. 
(@) and ( fl4|)) . Furthermore, the integrated weight of the 
function (pq) measures the deviation of the relaxation 
time th from the ordinary transport time Tt r - This de- 
viation is related directly to that of the Hall constant at 
zero frequency from its value at high frequencies: 



t h - r tr Rh{lo = 0) 



Ttr 



R 



(68) 



H 



In the phenomenological expression (|6g), the tempera- 
ture dependence is entirely contained in the two relax- 
ation rates. They are expected to vary as r fr cx 1/T 
and th oc 1/T 2 |J. For dimensional reasons, we take 
th oc J/T 2 , since J is the only energy scale character- 
istic of our model. From these properties, we may infer 
the following: For one thing, the memory function does 
not only describe the unusual frequency dependence of 
the Hall constant in high-T c superconductors. Also, the 
observed anomalous temperature dependence is mainly 
due to this memory function contribution. For another, 
we expect Tt r to be relatively smaller than rjj at low tem- 
peratures, since th/t^ oc J/T. Thus, Eq. ( |68| ) suggests 
that the Hall constant increases when zero frequency is 
approached. This enhancement was actually observed 
in the above mentioned measurements on YB^CuaOe+z 
p2fl . In any case, Eq. ( p8| ) indicates that the sign of the 
Hall constant is solely described by the high-frequency 
object R H , as claimed in Ref. ||. 

Does the discussion so far point towards spin-charge 
separation as advocated by Anderson? Obviously, the 
current operators J v are related to the charge degrees of 
freedom only. Consequently, the spin physics must be 
accounted for by the residual forces j v . This is also re- 
flected by the proportionality th oc J/T 2 . However, up 
to now, we do not have a compulsory argument why these 
residual forces should describe exclusively spin degrees of 
freedom. 

In summary, Anderson's notion of two distinct relax- 
ation rates is naturally backed up within the Mori theory. 
They may be interpreted as the time scales set by the 
current operators and their associated residual forces. 



B. Moments approach to the memory function 

Of course, it would be interesting to calculate the func- 
tion kiyi) of Eq. ( |64| ) quantitatively within the t— J model 
in order to relate its width and its integrated weight to 
the parameters t, J, temperature T and doping 6. From 
the above discussion, we expect that the relevant infor- 
mation about its overall form may be put into only few 
parameters. As already mentioned, even two parameters 



as in Eq. (66) have been sufficient to obtain an excel- 
lent fit of experimental data In this subsection, 
we suggest the following procedure to construct k(u): 
Parametrize this function by n parameters which are sub- 
sequently fixed by its first n moments. Up to now, this 
seems to be the only reliable way to take into account 
the superprojectors Q. Of course, all moments 



dto k(uj)oj 



21 



(69) 



exist, why we have to replace the Lorentzian ( p6[ ) by a 
"short-range" function, e.g. a Gaussian multiplied by a 
polynomial. We proceed by relating the moments of k(uj) 
to that of the function s(w) and the optical conductiv- 
ity. Finally, we discuss a possibility to calculate these 
moments. 

The ordinary conductivity is given in Mori theory by 
the following expression B: 



a xx {u + i0 + ) = ie 2 l3{J x - \ 



1 



uj + L + i0+ 



\Jx;o) 



(70) 



the real part of which may be shown to be an even func- 
tion of uj |^|, and, due to a sum rule, can be written in 
terms of a function c(lu), that is normalized to unity: 



U{a xx (uj + i0 + )} = 7re 2 x°c(Lu). 



(71) 



In the appendix, we show that all the moments ( |69| ) may 
be traced back recursively to that of the functions (p2 
and (fnl): 



si 



Q 



dui s{w)uj 21 
duj c{u>)u) 21 . 



(72) 
(73) 



provided I > X. Thereby, only the definition of the su- 
perprojectors Q has to be used in Eq. (|l7|). The result 
may be written as follows: 



si 



i 

3 = 1 



for I > 1 . 



(74) 



where the coefficients a,j are polynomials in the moments 
(|73|) and are listed in the appendix up to j — 6. If we 



had a good method for calculating the moments (72) and 
([73|), we could construct the unknown function k\uS) via 
its first n moments as already explained. This approx- 
imation is reliable, if the zeroth moment fco, calculated 
from Eq. (|6^), converges fast enough with increasing n. 
Since in the t — J model, k(uj) has only one peak around 
lu = 0, the first few moments are expected to be sufficient 
for this to happen. An advantage of this approach is that 
moments are global properties of a spectral function and 
as such are less sensitive to its detailed resonance struc- 
ture and to approximations involved. Also, some approx- 
imation schemes are better suited for the calculation of 
moments rather than the underlying spectral function. 
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For instance, within a high-temperature expansion, 
moments are accessible, at least in principle, while the 
corresponding spectral function is not. However, in the 
case of the moments (F^J), only very few moments, and 
to only low orders in l/T, are within reach. This is all 
the worse, since now, we are interested in the Hall con- 
stant at low frequencies, i.e. no other high-energy scale 
is present as it was in the context of R* H . 

Another example is the exact diagonalization tech- 
nique |p3[ . In this method, spectral functions are calcu- 
lated numerically via the exact eigenstates and eigenener- 
gies on the basis of their Lehmann representation. While 
being exact, the intrinsic problem of this method is the 
constraint of working on relatively small clusters. There- 
fore, the delta functions have to be broadened in or- 
der to obtain smooth functions. In contrast to this, 
no additional approximation is required when calculat- 
ing moments. Of course, the smallness of the clusters 
remains the major restriction of this method. Never- 
theless, when combined with the moments approach as 
suggested above, it is a means of extracting reliable in- 
formation about the frequency dependent Hall constant 
within the t— J model and should therefore be the subject 
of a future work. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have derived first results for the Hall 
effect in correlated electron systems in the strong correla- 
tion limit within the recently developed memory function 
approach |IJ|. We have focused our attention mainly on 
the memory function term, which is neglected in Boltz- 
mann type approaches. The important new physics to 
be expected from this contribution comprises the unusual 
frequency and temperature dependence of the Hall con- 
stant, observed in the normal state of high-T c supercon- 
ductors. 

In the single-band Hubbard model, a finite amount of 
spectral weight for particle-hole excitations, caused by 
the Hall current, is always pinned at the energy U. This 
is valid no matter how large the correlation strength U 
is. We have shown that the memory function removes 
these high-energy excitations in the limit U — > oo, thus 
accounting for the frequency dependence of the Hall con- 
stant down to frequencies within the charge transfer gap. 
The corresponding decrease of the Hall constant by a fac- 
tor (1 + 6)/2 was calculated exactly to leading order in 
l/T and corroborated within an approximate treatment, 
valid at T — 0. However, our analysis did not provide us 
with information about the frequency dependence of the 
Hall constant at lower frequencies. The reason is that 
it was based on moments, which are dominated by the 
high-frequency contributions. We have also calculated 
the infinite frequency Hall constant analytically within 
the so-called Hubbard I approximation. In essence, we 
recovered the exact result for U = oo and T — > oo and 



explained, why this result does not rely on the location 
of the Fermi surface. 

Finally, the Hall constant at low frequencies was in- 
vestigated within the t — J model, an effective model 
acting in the reduced Hilbert space without doubly oc- 
cupied sites. We observed that our memory function for- 
malism distinguishes inherently between two time scales: 
Firstly, the dynamics of the current operators is char- 
acterized by the ordinary transport relaxation time T tr - 
And secondly, the impact of all the other degrees of free- 
dom on this charge transport is taken into account by 
fluctuating forces, that introduce an unusual time scale 
th • On the other hand, it was pointed out by Anderson, 
that temperature dependences of transport and Hall ef- 
fect measurements can best be understood in terms of 
two relaxation times, following a l/T and l/T 2 law ||. 
We have shown that the time scales encountered within 
the Mori theory are identical to those introduced by An- 
derson. Furthermore, we have shown that the deviation 
of the unusual decay time th from the ordinary transport 
time T tr is intimately connected to the frequency depen- 
dence of the Hall constant at low frequencies. Thus, the 
temperature and frequency dependence of the Hall con- 
stant result from each other and are both due to the 
memory function contribution. It would be very inter- 
esting to investigate this interplay further, both theoret- 
ically and experimentally. As for the theoretical side, we 
have proposed an approach based on moments. It allows 
the exact treatment of the superprojector that reflects 
the distinction between the two time scales, leaving us 
with the problem of finding the first few moments of the 
ordinary current-current correlation functions. Except 
for a well studied prefactor, the memory function term 
is mainly determined by two parameters. Therefore, we 
expect its first few moments to provide us with enough 
information to fix these parameters. 
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APPENDIX A: REDUCTION OF THE 
MOMENTS OF THE MEMORY FUNCTION 

First of all, we recall the definition of the moments ( |l5| ) 
and (|l6|), however with all operators and superoperators 
now being redefined within the t — J model. They are 
related to the moments ([39]) and ([72]) via the equations 
I<i = -i < [J x -fl, J y ;a\ > h and Si = -i < \J X ;0, Jyfl] > 
s;. In addition, we require the moments ([73]), renormal- 
ized as Cj = x° c i- Next, we define the quantities 
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X? := i/3(L n J X ; \ LQ ...QL \J y , ) , (Al) 

v ' 

I operators L 
:=P(L n J X ;o\ LQ . . . QL \J X . ) , (A2) 
I operators L 

for I > 1. By using the definition (Q) for each of the 
first superprojector on the right hand sides, we obtain 
the following recursion relations: 



12c?c 3 



a 3 = 2c 3 - 6cic 2 + 4c? , 
(24 = 2C4 — 6C1C3 — + 12c? c 2 — 5c^ , 
05 = 2C5 — 6C1C4 — 6C2C3 

- 20c? c 2 + 6c? , 
ag = 2c6 — 6C1C5 — 6C2C4 

+ 24cic 2 c 3 - 20c?c 3 4 



12ciC2 



12c? c 4 - 3cj 



4rf 



30c? c^ 



30c?c 2 



7c? . 



(A19) 
(A20) 

(A21) 



(A22) 





vn+l 

~ A l-1 


1 1 v 

- ^-1 

- ^-1 


for n odd , 


(A3) 




XT 




for n even , 


(A4) 




yn 


_ v n+l 
— 


for n odd , 


(A5) 


[1] 




_ vn+l 

— 


for n even 




(A6) 


[2] 



To prove these relations, we only have to use the following 
facts: A™ and V" vanishes for all odd and even integers 
n, respectively. This is due to the fact that the functions 
( |l3| ) and (|7l]) are even, hence the ir odd moments van- 
ish. Moreover, the quantities ( |Al[ ) are of first order in 
the magnetic field. To relate the unknown moments Ki 
to the moments Si and Cu we have to supplement the 
recursion formulas (A3)-([Aq) by the following equations: 



K, 



Y° 

^21 + 1 1 



Si = -X? 



Ci=Y, 



21-1 



(A7) 
(A8) 
(A9) 



However, only even numbers of iterations occur. There- 
fore, we may combine two successive iteration steps into 
one. This leads to the following effective recursion rela- 
tions: 



C, 



S n 



"+ 1 ttO 1 u n atI 
L-0 L-0 



N n = N n+X_n N l_^ 
On 



where we have defined 



H? := X. 



2n 



N? 



20 + 1 
v 2n-l 
1 2l+l 



for n > 0, I > , 
for n > 1, I > 0, 



(A10) 
(All) 



(A12) 
(A13) 



and where the contact to the moments Ki, Si and Ci is 
established by means of the equations 



Ki 

Si = -Hq , 



Q = Nn 



This recursive procedure results in Eq. ([? 
six coefficients being given as follows: 



ai = 2ci , 
a 2 = 2c 2 - 



3c?, 



(A14) 
(A15) 
(A16) 

with the first 

(A17) 
(A18) 



[5 

[e; 

[7; 

[8. 
[9 

[10 

[11 
[12 
[13 
[14 

[15; 



[16 
[17 
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[19 
[20 

[21 

[22 
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